УДК 550.931--550.838 
Відновлення контактної границі в шаруватому середовищі 


О Ю.І. Дубовенко, 2002 


Інститут геофізики НАН України, Київ, Україна 
Надійшла" 7 червня 2002 р. 
Представлено членом редколегії В.І. Старостенком 


На класі Нумерова нелінійну контактну задачу з задовільною точністю апроксимовано розв'яз- 
ками лінійних інтегральних рівнянь першого роду двох типів. На їх основі зформульовано задачі ко- 
ректного відновлення контактної поверхні за значеннями поля на майже необмежених і істотно об- 


межених множинах. Для розв'язання цих рівнянь запропоновано стійкі ітераційні процеси та схему 
прискорення їх збіжності. 
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При інтерпретації даних потенціальних полів з метою вивчення глибинної будови Землі 
визначальною ланкою є розв'язок зворотніх структурних задач потенціалу | 1-4). Шуканими 
параметрами в них є характер поведінки й глибина границь розділу шарів, на яких функції 
щільності (намагніченості) зазнають розриву. Математично розв'язок таких задач зводиться 
до розв'язку інтегральних рівнянь 1-го роду й у цілому нестійкий. 

Загалом нелінійну зворотню задачу логарифмічного потенціалу для контактної поверх- 
ні зведено до розв'язання нелінійного інтегрального рівняння типу Урисона 


гі 2 2 
ПЕ) мод () 
ло (5-х) 1075) 2л/з 
яке відрізняється від рівняння Б.В. Нумерова |3| наявністю замість середньої глибини й до 
контакту значення 0 (х) цього контакту у фіксованій точці х; про можливість такої підстанов- 
ки згадував О.А. Шванк |б6). Розв'язок наведеного рівняння знайдено на компактному струк- 
турному класі - так званому класі Страхова - перетині класу функцій щільності, обмежених 
і локально інтегрованих майже всюди в необмеженій замкнутій області, що містить тяжі- 
ючі маси, і класу неперервних за Гельдером (тобто обмежених) функцій контакту; на цьому 
класі розв'язок нелінійної контактної задачі (1) існує, єдиний 1 стійкий |3, 4). Поведінку рів- 
няння (1) в комплексній площині вивчено в працях |7, 3). Одначе, в задачах типу (1) точний 
розв'язок досягається винятково при заданні значень поля на необмежених інтервалах, тоді 
як у польових умовах елементи поля вимірюють на істотно коротких профілях. Отже, актуа- 
льним є вивчення контактної задачі за умов задання вхідних даних на істотно обмежених 
множинах 
В зв'язку з цим розглянемо лінеаризований варіант задачі |9), який дозволяє віднови- 
ти контакт за значеннями поля на істотно обмежених інтервалах за допомогою заміни інте- 
грального рівняння з повільно спадаючим ядром Пуассона рівнянням зі швидко спадаючим 
ядром. Наведемо способи її коректного розв'язання з врахуванням результатів |7, 10, 11). 
Зазначимо, що поняття "істотно обмеженого інтервалу" запозичено з роботи |12) з 
наступним уточненням: нехай потрібно визначити в рівномірній дискретній мережі точок х, 


з постійним кроком Дх деякий фрагмент с(х)х є Їсуьс»| контактної поверхні сх) за виміря- 
ними з відомою точністю ІнС2)- ца (о «5 на інтервалі Где, ТЕ сьо | значеннями поля и(х). 


В процесі розв'язання цієї задачі з допомогою деякого ітераційного процесу |3| на основі 
рівняння (1) потрібно обчислити таке значення межі інтегрування а прямого оператора (1), 
при якому похибка його апроксимації не перевищить величини, кратної похибці поля 


Іо) єс а)| « 15, 125, де «(х, а) - точність скінченновимірного проектування, «(х,п) не 


точність чисельного інтегрування (додатковий член формули чисельного інтегрування). Як- 


що при цьому поле и, (х) задано на інтервалі Гор, с так, що виконується вкладання 
Ї сь офедх 2,6, 3ач- Й! 2| 2 со, то інтервал будемо називати майже необмеженим, в ін- 
шому разі - істотно обмеженим. У першому випадку за співвідношенням (1) непогано від- 
новлюється контакт, у іншому - надійно відновити контакт з допомогою ітерацій з ядром 
Пуассона неможливо і необхідно застосовувати інші, більш швидкозбіжні схеми обчислень!. 
Моделлю геологічної будови для вищезгаданої задачі є два зоряні відносно нескінчен- 
но віддаленої точки однорідні циліндричні шари різної щільності с, що залягають один над 
одним й обмежені смугою 
Пеїкхбьоо«х«оу, «С(х)є А), (2) 
де зіпі сх) , П" з8ир сх) - глобальні екстремуми контакту 4Є37 а постійні й, - ші г) і 
йо зару | описують множину контактів с(х) зі смуги (2) за умови 0 « Й, 5 сх) «ль ХЕ Ек 
Якщо в класі Страхова 5, п) І9| виділити на підставі співвідношень 
ос) со), 2) - й ух є В, 3) 
0.6п «(х) «141, О« ПР «АЙ", (4) 


що обмежують "коливання" 71 Га ) контакту, деякий компактний клас Нумерова Ми (а, п) 


то досить легко показати, що для контактних границь, описаних функціями з цього класу, 
головну роль у гравітаційному притяганні підземних мас відіграє перша складова в лівій 
частині рівняння (1), а друга - лише роль деякої малої поправки до неї. |11|, Наслідком цьо- 
го твердження є можливість подати на класі Ми (1, п) нелінійне рівняння (1) з точністю 


22 (х) б ле) - "Р лінійним інтегральним рівнянням Фредгольма 1-го роду з ядром Пуассона 


| ян (244 а (хо) ех) Лух є А, б 


причому розв'язок Сх), х є ВО цього рівняння на класі Нумерова існує і єдиний при заданні 


|» о) 


правої частини на довгому профілі. 

Досліджувана задача полягає у відновленні контакту при обмеженнях (2-4) за полем, 
заданим на "необмежених" |с,, с; | і істотно обмежених |с|ьс, | множинах|. Для відновлення 
контакту при заданні поля на істотно обмежених множинах рекомендується еквівалентне 
рівнянню (5) лінеаризоване рівняння зі швидко спадаючим ядром Шварца (111: 





5 б пе бує -мх), (6) 
зануа | (а мезо а, пуб 


Його розв'язок існує (як наслідок лінеаризації) і єдиний. Однак це не єдиний шлях визначен- 
р У р ц Д 
ня наближень 0 (А): при належності шуканого розв'язку до компакту 0 (х)є Ми") (1, п) лі- 


нійне наближення контакту 0 (2, п) отримуємо як суму розв'язків 0 (п) зі (ж, п) 5 "(ол) 
лінійних інтегральних рівнянь 


зна маних кі стае) 


Існує ще один варіант знаходження наближення : (2,1) як розв'язку рівняння 


Коло 28 (ол) | Гео Об (8) 


1 Одну з таких схем для рівняння (1) якраз і наведено в роботі (12). 


З 


де 5 "051) визначено з першого рівняння системи (7). Такий альтернативний розв'язок від- 
різняється від точного розв'язку 0 (х) рівняння (1) не більше за квадрат відхилення 
С (х) б лог) - А, «АХ ,хє ЕФ), однак він значно гірший за якістю, ніж наближення, от- 
римані за формулою (7) через наявність інтегралу з повільно спадаючим ядром Пуассона. 
Розв'язки рівнянь (5-8) залежать від апріорі невідомої середньої глибини до контакту ЛП (що, 
до речі, є наслідком модельних зображень складних середовищ, а не принциповою ха- 
рактеристикою задачі). Її значення з геологічних міркувань не завжди відомі надійно, тому 
важливо розробити спосіб знаходження контакту, у якому було б усунуто подібну парамет- 
ричну залежність. Цього можна досягнути шляхом введення нових узагальнених співвідно- 
шень для контакту 





19 ру і -и(х), сх) є (4)- 6), (9) 
НО даху 9 


23 сз 
де моно У о му, 
пт 0 3 (2К0(х)) (5-х) 
Подібна заміна справедлива й для рівнянь (7) 1 (8). Практична значимість цих виразів полягає 
в можливості отримувати значення глибин контакту в кожній фіксованій точці рівномірної 
мережі в ході ітераційного поліпшення моделі, що надійніше, ніж підстановка апріорних зна- 
чень й. Співвідношення (5-6, 9-10) цілком працездатні й у випадку магнітоактивних контак- 
тів, коли в їхніх правих частинах стоятимуть похідні напруженості магнітного потенціалу; 
контакт відновлюють із системи спряжених рівнянь типу (7) як суму значень 

цу) є 57 (х, п) 57 (о А) м, А) 57 (х, п) 57 (ол). 
Ця обставина відкриває перспективи для комплексної інтерпретації гравімагнітних полів і 
вимагає ретельнішого вивчення. 

Питання розв'язності рівняння (5) досліджував Б.А. Андрєєв. Для його розв'язання 
він запропонував спосіб послідовних наближень (13, а М.М. Лаврентьєв (14) обгрунтував 
його застосовність для розв'язку операторних рівнянь першого роду?. Згодом розв'язність 
рівнянь (3) 1 (9) вивчалась в роботі |3|, в якій для їх розв'язання було запропоновано так звані 
ітераційні процеси Лаврентьєва- Андрєєва 


мне у Ен РО)оМо, 0 аю 





1 с (1) Х 
ОР ОРУ О  УРРУР А МЛ РЕ 
я (СМ) (6-х) 
Розв'язність рівнянь (6, 10) вперше дослідив А.В. Чорний |3|; для обчислення їх розв'язків 
згодом введено нові ітераційні процеси, названі процесами Лаврентьєва-Чорного 


«уч еВ арени) РО) б здано 0) 





а 





оче бо СОУеаЕ чех), ех) є м) є 0,2,со | (14) 





2 Ці обставини відбито у назві наведених нижче ітераційних процесів. 
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Дослідження збіжності цих процесів здійснюється елементарними засобами (методом 
індукції за умови збіжності мажоруючих часткових сум) внаслідок еквівалентного подання 
розв'язків через наступні наближення: 


СО) а їмеме| 2 (І -е"й у сова - ха, (15) 
2 Імені У (і я (освожу сова - хо, па 0, со (16) 
о 0 ко0 


Твердження про збіжність цих наближень водночас, згідно пропозиції (15|, є конструктивни- 
ми теоремами існування розв'язків рівнянь (3) 1 (9); аналогічні співвідношення, лише склад- 
ніші в плані чисельної реалізації, виведено для узагальнених рівнянь (6) 1 (10). Вивчення пи- 
тання про швидкість збіжності різних ітераційних алгоритмів? для рівнянь (15) 1 (16) дозво- 
ляє дійти висновку, що послідовні наближення Лаврентьєва-Андрєєва у ГЕ); і Лаврентьєва- 
Чорного у ЕТ рівномірно збігаються до граничної функції 0 (х) є Ми) (1, п), яка задові- 


льняє кожному з двох еквівалентних рівнянь (35) 1 (9) зі швидкістю геометричної прогресії, 
оцінюваної , відповідно, величинами 


М ео ки ЗА 


Р ОР ПР Че ЛИ 
пацдеЗ би (хх На-к2)! Нана); 
а модулі неперервності зворотніх операторів цих рівнянь оцінюються як 


кчено є н косе), 
де М (п) щ 2тах|(х); м) щ 2 тах|у(х,21); МУ) ш 2 тах|мх) є мх21) , хе А а м(х, 21) 


і мо2л) - аналітично продовжені значення м(х) і мох) як гармонійних функцій з рівня у - 0 

















на рівень у - 2п. Схожий результат справедливий і для узагальнених рівнянь контакту. 

Порівняння модулів неперервності для рівнянь (5) і (6) демонструє зворотну пропорцій- 
ність швидкості збіжності ітерацій до глибини Й і майже однаковий її порядок для обидвох 
процесів за невеликої (усього вдвічі) переваги у швидкості збіжності процесу (13) над проце- 
сом (11). В обох процесах швидкість спадає з наростанням глибини Й й індексу ітерації п, але 
практично ітераційний процес Лаврентьєва-Чорного має перевагу над "андрєєвським": завдя- 
ки швидкому згасанню по горизонталі ядер перетворень (13) для відновлення наближень 
5(х, А) з визначеною точністю можна обмежитись значеннями поля и(х) на істотно корот- 
ших профілях, ніж при використанні будь-яких схем з ядром Пуассона. Одначе, останніми за 
певних умов можна користуватись, якщо істотно прискорити їхню збіжність. 

У процесі вивчення питань можливості розв'язання рівнянь (3) і (9) попутно було по- 
ліпшено формули, що уточнюють наближення Нумерова (1), отримані колись Б.А. Андрєє- 
вим (13), М.Р. Малкіним, А.К. Сеньком 1 А.К. Маловичком (16-18). Реального поліпшення 
ітерацій Б.А. Андрєєва можна домогтися, якщо лінеаризоване наближення розв'язку неліній- 


ного інтегрального рівняння для контакту 0 (х)є Ми (на (1, п) апроксимувати виразом 


п-1 3 
бо вано) У (-уночя |оР 0) - (20) 17 
5 (, ) ц(х) А! ) па2 Гаучочх) Їм(є ) цх ) 5 ( ) 
який від істинного ухиляється не більш, ніж на величину лоЧЄ) при Л-»о. 


Широко вживані в практиці аналітичного продовження формули Нумерова, Хьюза, Стра- 
хова є частковими випадками формули (17). Дійсно, якщо брати за нульове наближен- 
ня вираз 


З Відсутність таких даних заважала оцінці їх чисельної ефективності і можливості порівняння різних алгоритмів з метою 
вибору оптимального для розв'язуваної задачі. 
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о 


1 1 П 
х)а ---чи| х,0)- - | ---онноно ЩО, Й) - м х, А/С 2, 18 
Р ав 
де п - відстань до найближчої особливості поля и(х,0), неважко отримати для індексів 
п - 0,1,2 наближення Б.В. Нумерова | 31, Д. Х'юза (181, та В.М. Страхова (19), відповідно, 


пен ро у МО, Ах) 2), 


ОР З ВО У О.Б У 
Сол) (со) у ий (ха, (19) 


02) М 17 Й 2Й ЗИ 
зо. ї гору с М Су С неуеої зо 
Ядра цих перетворень додатні, а ядро для п - 4 - уже ні, що призводить до осциляції 
розв'язку, тому в практичних обчисленнях варто обмежитися для знаходження контакту од- 
нією з наступних наближених форму 
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КОТРА У ВАРИ РУ -! а в бо з м дид Є, ди зи(2)- ц(х) 
націо ПЗВ ІЗ В хо? (6-х) , 
інно ду Йайчє Рів нстічй 0 09 
-х) «2917(5 - х) «(Е-х)"|й 












Диас. 


тт 


зн) (вх) ак) (6-х) 

Суть пропозицій М.Р. Малкіна й А. К. Сенька стосовно до постановки задачі (2-5), які 
уточнюють розв'язок рівняння типу (35), зводиться до наступного. Малкін запропонував уточ- 
нювати кожне з лінеаризованих наближень (20) контакту 0 о, п) за наступною схемою 


Со 1) з С л)ч па ря - в якій ми обу за величину поправки брати 


пе о г 5 (ел) євьнурав,. 


Ця поправка підвищує точність кожного з зо В наближень на порядок щодо початко- 





вих - аж до величини й (п). А.К. Сенько для уточнення кожного з лінеаризованих набли- 


жень (20) розв'язку 0 еп) нелінійного рівняння для контакту (1) запропонував наступну 





послідовну схему 2" (х, )- Я а ВТР АС), а 0,1,2.3,...; Якщо поправку Д/(х) знахо- 


дити як розв'язок наступного лінійного бо рівняння 1-го роду 
і уча 
пАС" (1 6) 


13. (е-ху я ІСОСЮЇ 
(нд | еоічрння ВЕ СР 

ж (вч) НО) 
то ітерації будуть генерувати наближення, що ухиляються від точного розв'язку не більше, 
ніж на величину 05 / (р) бот тах| АС (х). А.К. Маловичко |18| запропонував природнє 
узагальнення ітерацій Б.В. Нумерова й А.К. Сенька для уточнення контакту, яке на підставі 
ітераційного процесу (11) нам вдалося удосконалити. Пропонується новий спосіб послідов- 
них уточнень контакту за схемою 2, (х, й) є 2, (х, п) АС, (ж, А) «А «А, у вигляді процесу 


Аа) ці | пу авто) 01) 





с 2 2 
де АС) - сО)- 6 (о п) - Па Ізеоі но) 48, тео, пеб,со, 
жі (5-х) кСЦ5Л) 
ох) я Сх В) о 01,23. Ал. 
За умови 0 (х)є Ми и П), х є В? послідовні наближення, генеровані процесом ітерацій 
(21), збігаються при п -» оо до граничної функції 0 (х) зі швидкістю геометричної прогресії, 
отже, існують. 

Висновок. Точність обчислень контакту залежить не лише від похибок поля 96, але їі 
від довжини інтервалу його задання; за однакової точності обчислень у випадку задання 
значень у(х) на короткому інтервалі для відновлення їх за допомогою операторів зі швидко 
спадаючими ядрами потрібний менший обсяг інформації, ніж у випадку задання поля на не- 
скінченних профілях. Строге обгрунтування запропонованих методів відновлення контакту 
при обмеженому профілі вимірів дозволяє рекомендувати їх для автоматизованої інтерпрета- 
ції даних гравіметрії в маловивчених районах. 

Виношу глибоку подяку В.І. Старостенку і А.В. Чорному за критичні зауваження, які 
посприяли поліпшенню й глибшому розумінню матеріалу цієї статті. 
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